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Abstract 

The resummation of 0{a l ^ rl \r? 1 x) terms in the evolution equation of the singlet part 
of g±(x,Q 2 ) is carried out. The corresponding singlet evolution kernels are calculated 
explicitely. The leading small- x contribution to the three-loop splitting function matrix is 
determined in the MS scheme. Relations are derived for the case of N = 1 supersymmetric 
Yang-Mills field theory. Numerical results are presented for the polarized singlet and 
gluon densities, and the structure functions gf{x,Q 2 ) and gi(x,Q 2 ). They are compared 
for different assumptions on the non-perturbative input distributions, and the stability of 
the results against presently unknown subleading contributions is investigated. 
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1 Introduction 



The small-x behaviour of polarized structure functions is a largely unexplored subject. The 
current measurements cover at most the range x > 5 • 10~ 3 |TJ. At even lower x- values both 
size and sign of the structure functions gf(x,Q 2 ) and g™(x,Q 2 ) are yet unknown. The initial 
distributions at a starting scale Qq of the QCD evolution can not be determined within perturba- 
tive QCD. Irrespectively of the specific behaviour of these non-perturbative quantities one may 
consider, however, the QCD evolution of the structure functions and parton densities. Besides 
the well-known leading order anomalous dimensions ||, recently also the next-to-leading order 
ones have been calculated . Furthermore the all-order resummation of leading singular terms 
in higher-order anomalous dimensions may be important in the small-x range. These terms 
behave as a 1 ^ 1 In x, corresponding to N(a s /N 2 ) l+1 for iV — ► in the Mellin moment plane. 
Their numerical effect on the evolution of non-singlet [[| structure functions has been analyzed 
in refs. |5|, ||. Contrary to an earlier expectation in ref. |7j it turns out to be very small and is 
found to depend on the way in which fermion number conservation is imposed, i.e. also on less 
singular In x contributions. Recently also an equation for the leading small-x resummed singlet 
evolution has been obtained 0. 

In the present paper we study this singlet resummation in the framework of the renormaliza- 
tion group equation. After recalling the general evolution equation and setting up our notations 
in Section 2, we calculate the anomalous dimension matrix as a series in a s explicitely in Sec- 
tion 3. A series of properties of this matrix is discussed. The numerical effect of the resummation 
beyond the known |J next-to-leading order effects on the small- x behaviour of the polarized 
singlet quark and gluon densities, AT>(x,Q 2 ) and Ag(x,Q 2 ), and on the structure functions 
gl (x, Q 2 ) and gi(x, Q 2 ) is then studied in Section 4, using input distributions as determined in 
recent analyses . We also investigate the stability of the results against yet unknown effects 
of terms less singular as x —> 0. Section 5 contains our conclusions. 



2 The Evolution Equation 

The evolution equation for the polarized singlet quark and gluon densities (AS, Ag) is given by 



ainQ2 \ Ag(x,Q 2 ) ) v->"»<~ y Ag(x,Q 2 

Here <8> stands for the Mellin convolution, and the matrix of the polarized singlet splitting 
functions P(x, a s ) is specified below. In the following, we will simplify the notation by using 
the abbreviation a s = a s (Q 2 )/47r for the running QCD coupling. The scale dependence of a s is 
obtained from 

-f>" +2 A, (2) 

fc=0 



d\nQ 2 



where only (3 = (11/3) C A - (4/3) T F N f and ft = (34/3) C\ - (20/3) C A T F N f -AC F T F N } enter 
up to next-to-leading order (NLO). Here C A = N c = 3, C F = (N 2 - 1)/(2JV C ) = 4/3, T F = 1/2, 
and Nf denotes the number of flavours. The matrix P(x, a s ) can be represented by the series 

oo 

P(x,a s )=J2a l s +1 P {l) (x). (3) 

1=0 



1 



Unlike the LO and NLO (MS) transition matrices P (0 \x) and P (1 \x) §, g], the splitting 
functions P^ >2 ) (x) are not completely known so far. For the solution of eq. ([3]) beyond NLO in 
Section 4, we will use the asymptotic small-x form of the latter matrices as derived in Section 3. 
From that solution the structure function gi(x, Q 2 ) is finally obtained by the convolution 



p,n I 

9i 1 



(x, Q 2 ) = ± {c 1>g {x, Q 2 ) AE(x, Q 2 ) + c 1>9 {x, Q 2 ) ® Ag{x, Q 2 )} + g^ s (x, Q 2 ) (4) 
with 

AS = An + An + Ad + AS + As + As , (5) 

and Am, Ad, As denoting the polarized up-, down-, and strange-quark distributions. We con- 
sider only the contribution of the three light flavours. The non-singlet part of gi(x, Q 2 ) has been 
dealt with in refs. ||, |6| already, to which we refer for further details. The coefficient functions 
Ci(x, Q 2 ) can be expanded in the strong coupling as 

oo 

c M (x, Q 2 ) = 6(1 - x)5 iq + £ aicg(x) . (6) 

i=i 

Here it is important to note that in the MS scheme the known coefficient functions c[ (x) for 
both I = 1 and I = 2 (cf. refs. 0, PH ) behave only like 

c?oca s ln(-) , ^ oca 2 s \n 3 (-) (7) 



X J \x 



at small x. Therefore a prediction can be made on the small-x behaviour of the three-loop tran- 
sition matrix in the MS-scheme, P^Lqi see ec L- Ql4 ) below. Note that in eq. (P the resummation 
under consideration is of leading order, i.e. the respective terms emerge only together with /3 
after rewriting the evolution equation in terms of a s . Details of the solution will be presented 



elsewhere 14 



3 Resummation of dominant terms for x — > 

The resummed transition matrix for the leading singular terms as x — > 0, P(x, a s ) x ^ , can be 
obtained from the solution of eq. (§) via inverse Mellin transformation: 

OO 1 

P{x,a s ) x ^ = £ P^l a l s +1 ln 2l x= ^M' 1 [F (N,a s )} (x) . (8) 
i=t 8tt 2 - - 

The matrix-valued function F (N,a s ) is subject to the relation 

F (N, a s ) = 16vr 2 ^ M - ^ F S (N, a s ) G + JLi F 2 (N, a s ) (9) 
derived in ref. ||, where Fs(N,a s ) is the solution of 

F 8 (N, a s ) = 16tt 2 ^ M s + ^C G -±- F 8 (N, a.) + ~ F 2 S (N } a.) . (10) 

The basic colour factor matrices are given by 

M ( C F -2T F N f \ ( C F \ (C F -C A /2 —T F Nf \ 

Mo = [2C F AC A J ' Go = C J ' Ms = C A 2C A ■ 



We determine the matrix P(x, a s ) x ^o in terms of a series in a s since this representation is needed 
for the solution of the singlet evolution equation ([]]), see Section 4. 
The lowest order entries read 



,(0) 



C F -2 T F N f 
2C F ACa' 



:i2) 



>( i) ( C F {2C A -3C F -4T F N f ) -2T F N f (2C A + C F ) 
^ 2C F (2C A + C F ) 8 C\ - 4 T F N f C F 

21 



(13) 



They agree with the respective leading In x terms of the complete LO and NLO splitting 
function matrices [@, ||. We also list the entries of P^o an d P^o explicitely in the colour 
factors: 



-5 C% - - C\ + 6 C A C F - 8 T F N f C F -Q T F N f C A 



P(2) _ 

99 


Ho 

3 


p(2) _ 
99 




p(2) _ 
99 




P(2) _ 
99 


2 p 

3 " 



'J Uf — VI F l\f 

z, 

, ,\ , -15 C\ + 2 C 2 F - 6 C F C A + 8 T F N f C F 

15 - 2 Cf + 6 C F C A - 8 T F N f C F 
28 + 2 T F iV/ - 4 T F Nf C\ - 24 C F T F Nj C A 



(14) 



and 



P(3) 
99 



p(3) 
99 



p(3) 
99 



p(3) 
99 



^o F 

45 



1 Q 

. 4 - 20 C F C A + 22C A C F -—C F - 74 T F N, Cj - 44 C F T F Nj C A 



'2 i /in /T~> at *\2 



—T F N f 
45 ' 



+ 2 TpiV7 C F + 40 {TpN f) C p 
-54 CV - 2 C 2 F C A + 40 TpNjCp - 128 - 8 T F N f C A + 15 Cf 
+ l08T F N f C F CA 



45 



2 

45 



54 C F Ci + 2 C F C A - 40 T F N f Cp + 128 C\ + 8 T F N f C A - 15 C F 
- VQ%T F N f C F C A 
-288 T F Nf C F C\ - 64 CpT F N f C A + Q T F N f C 3 F + 40 (T F N f ) 2 C F 

f w 



20 T F iV f Ci + 252 Ci 



(15) 



Comparing eqs. (|T2 HT5|) it is interesting to note that the leading small-x off-diagonal elements 
are related by 

PS/(T F N f ) = -Pg/C P . (16) 

We have verified this property analytically up to order I = 100. f] The numerical values of the 
matrix elements of jP^I+o up to Z = 10 are listed in Table 1 for SU(N C = 3), leaving the number of 
(massless) partonic flavours as the only free parameter since the corresponding analytical results 
become rather lengthy. The matrix elements for even higher indices are easily obtainable. 

2 We used the program system Maple V (lj] for checks and the derivation of higher order coefficients given 
subsequently. 
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N f 


= 3 




1 


P (0 


P (0 

r Q9 


P (0 
r g<i 


P (0 

r 99 




1 

2 
3 
4 
5 
6 
7 
8 
9 
10 


0.2666666667D1 
-0.1066666667D2 
-0.3679012346D2 
-0.6642085048D2 
-0.6110486315D2 
-0.3858083350D2 
-0.1679581048D2 
-0.5631967112D1 
-0.1424435573D1 
-0.2991318204D0 
-0.4868787168D-1 


-6.D0 
-44.D0 

-0.1394444444D3 
-0.2288296296D3 
-0.2154469209D3 
-0.1347153415D3 
-0.5955393524D2 
-0.1979044200D2 
-0.5081773820D1 
-0.1049686114D1 
-0.1756718578D0 


0.5333333333D1 
0.3911111111D2 
0.1239506173D3 
0.2034041152D3 
0.1915083742D3 
0.1197469702D3 
0.5293683133D2 
0.1759150400D2 
0.4517132284D1 
0.9330543234D0 
0.1561527625D0 


24.D0 
128.D0 

0.4188888889D3 
0.6981037037D3 
0.6685486038D3 
0.4201415122D3 
0.1868497637D3 
0.6213345892D2 
0.1601700122D2 
0.3302769905D1 
0.5557315074D0 






N f 


= 4 




I 


P (0 


P (0 

r Q9 


P (0 

r 91 


P (0 

r 99 



1 
2 
3 
4 
5 
6 
7 
8 
9 
10 


0.2666666667D1 
-16.D0 

-0.4953086420D2 
-0.8573278464D2 
-0.7633439480D2 
-0.4648843309D2 
-0.1955536907D2 
-0.6326464875D1 
-0.1545655016D1 
-0.3133281583D0 
-0.4922813282D-1 


-8.D0 

-0.5866666667D2 
-0.1788148148D3 
-0.2859456790D3 
-0.2595199393D3 
-0.1568583694D3 
-0.6687789759D2 
-0.2148071223D2 
-0.5319199339D1 
-0.1062889464D1 
-0.1712782960D0 


0.5333333333D1 
0.3911111111D2 
0.1192098765D3 
0.1906304527D3 
0.1730132928D3 
0.1045722463D3 
0.4458526506D2 
0.1432047482D2 
0.3546132893D1 
0.7085929761D0 
0.1141855307D0 


24.D0 

0.1226666667D3 
0.3905185185D3 
0.6315654321D3 
0.5831456986D3 
0.3540380256D3 
0.1519476618D3 
0.4881870016D2 
0.1214527167D2 
0.2420644158D1 
0.3928177051D0 



Table 1: The elements of the coefficient matrices in eq. (||) for Nf = 3 and Nf = 4. 

We have also calculated the matrices P^o f° r an A/" = 1 supersymmetric Yang-Mills field 
theory, i.e. Ca = Cf = 1, Nf = l,Tp = 1/2. One finds that in this case the so-called 
supersymmetric relation 

P%{x) + P$!(x) ~ PgW - P£(x) = (17) 

is fulfilled for the small- a; leading terms. We have verified this behaviour explicitely up to order 
I = 100. Beginning with 0{a 2 s ) even 

holds, and the matrix of the small-x transition functions depends only on one single scalar 
coefficient pi at each order in a s . Hence one can write 

oo 

P^sy = 2 fls Mi + a [ +1 ln 2/ x Pl M 2 (19) 
i=i 

with 

M^Mr Y ={\ ~\), M 2 =(-\ "M. (20) 



Note that there is a series of relations between the matrices Mi and M 2 : 



M 1 M 2 = 3M 2 , M\= M 2 , [Mi, M 2 ]=0, (21) 

which determines the structure of (|H5p, and a similar result is obtained for Ff VSY (N) . 

We have furthermore investigated whether the dominant terms found in the present smail-x 
expansion can be obtained in the large-Afy expansion [13) partly, where predictions have been 
obtained on the behaviour of the eigenvalue 

e 2 = \ [P„(N) + P gg (N)} - \sI\pJn) - P 99 {N)f + ±P q9 {N)P aq {N) (22) 

of the Mellin transformed all-order singlet transition matrix. However, as it is found for the 
lowest order in a s already, the contributions to the large Nf expansion of e 2 can be obtained 
only by accounting for non-leading terms for x —>■ 0. Therefore the limits x — > and Nf — > oo 
do not interchange and a further test on the elements of P^ l \N) can not be obtained in this 
way0. 



4 Numerical results 



After transformation to Mellin moments, the singlet evolution equation ([!]) is reduced to a 
system of coupled ordinary differential equations. Unlike in the non-singlet case considered in 
refs. || [| the solution cannot be given in a closed analytical form beyond LO here, due to the 
non-commutativity of the matrices P^\N) for different orders in a s . Instead the evolution 
taking into account the leading small-x resummed kernels (|8|) has to be written down in terms 
of a power series in a s , yielding 



AE(N, a s 
Ag(N, a s 



l + 5>iU«(JV) 



a 



i + £4e/«(at) 



i -i 



AS(iV,a ) 
Ag{N, oq) 



(23) 

with ao = a s (Qo). The singlet evolution matrices U^(N) can be expressed in terms of the 
splitting function moments P^^^N). Technical details can be found elsewhere 0]. Due to 
the structure of eq. (§), the solution ( p3|) is necessarily related to an asymptotic expansion. 
For all practical cases, say x > 10~ 6 , however, retaining 8-10 terms in eq. (|23| ) is adequate 
for obtaining accurate and stable results. The transformation of the outcome back to x-space 
finally affords one standard numerical integration in the complex iV-plane. 

We study the numerical consequences of the resummation using the recent NLO 
parametrization of ref. (GRSV) as the input at the reference scale Ql = 4 GeV 2 . The presently 
large ambiguities due to the virtually unknown polarized gluon distribution will be briefly illus- 
trated by the maximal and minimal Ag scenario of the same group. For a short review on current 
parametrizations of polarized parton densities see, e.g 
the /3-function is fixed at Nf 



ref. 15 



The number of flavours Nf in 
4 for all results shown below, and A = A^iV/ = 4) = 200 MeV 



is employed in the standard approximation to the NLO running of a s 



a s (Q 2 ) 



1 



/3 ln(Q 2 /A 2 ) 



Pi lnln(Q 2 /A 2 ) 
Pi ln(Q 2 /A 2 ) 



(24) 



^Recall that the non-singlet transition functions in the small- a; limit do not depend on Nf at all j5j. 
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The number of partonic flavours in the splitting functions (N) is restricted to Nj = 3 |16 . 

We are now ready to present the resummation effects on the polarized singlet quark and gluon 
densities, xAS and xAg, as well as on the proton and neutron structure functions, g[ and g™, 
respectively. The structure functions have been obtained from the (scheme-dependent) parton 
distributions by the convolution (|4]), i.e. without subtracting the resulting subleading a 2 terms. 
The non-singlet case has been studied for our present input choice in ref . || . The corresponding 
resummation effect is negligible (about 1% or less) over the full kinematical region considered 
here. In Figure 1 the NLO and leading small- x resummed (Lx) results are displayed for the 
singlet parton densities, and Figure 2 depicts the corresponding proton and neutron structure 
functions gf ,n . 

The resummation effects are much larger than for the non-singlet quantities, as to be expected 
from the comparison of the coefficients in Table 1 to the corresponding non-singlet results. E.g., 
the ratio (NLO + Lx)/(NLO) amounts to about 1.72 (1.64) for AS (Ac/), respectively, at Q 2 = 
10 GeV 2 and x = 10 -4 . It should be noted in this context that the small-x evolution strongly 
depends on the practically unknown gluon input distribution. This is illustrated in Table 2 where 
the resummed results of Figure 1 are compared at two representative values of x and Q 2 to those 
obtained by evolving in the same way the 'minimal Ag' and 'maximal Ag' distributions of ref. 0. 
One finds variations up to a factor of almost 5 (10) for AS (Ag), respectively, indicating that 
at present the input ambiguities are the dominant source of uncertainties also at small x. 



Q 2 


10 GeV 2 


100 GeV 2 


X 


10~ 4 


10~ 3 


10" 4 


10~ 3 


xAS 


-0.0100 
-0.0285 
-0.0473 


-0.0169 
-0.0396 
-0.0560 


-0.0171 
-0.0505 
-0.0855 


-0.0218 
-0.0523 
-0.0772 


xAg 


0.019 
0.101 
0.201 


0.034 
0.152 
0.294 


0.053 
0.226 
0.432 


0.071 
0.281 
0.528 



Table 2: A comparison of the resummed evolution of the polarized parton distributions for different 
assumptions on the gluon distribution Ag. Upper lines: minimal gluon, middle lines: standard set, 
lower lines: maximal gluon (and corresponding quark distributions) of ref. || at Ql = 4 GeV 2 . 

An obvious question concerning the large effects found is whether the resummed a 1 ^ 1 In 2 ' x 
terms really dominate with respect to the presently yet uncalculated terms less singular in lnx. 
Recall that we found for the non-singlet structure functions that this is not the case || 
Lacking any further information on the higher-order splitting functions, e.g. from sum rules 
as in the non-singlet and unpolarized singlet cases, it appears reasonable to assume that the 
coefficients of the first less singular term is of roughly the same size but of opposite sign as the 
leading one. To obtain a first estimate we use 

AP(' J>1 ) -» AP (i>1) • (1 - N) . (25) 

This assumption is motivated by corresponding relations in the LO and NLO splitting functions 
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where, e.g., for the largest quantity P gg : 

p£Hm = | - 15 + <w > p^w = 0(1/Ar) (26) 

for three flavours. The corresponding results are also given in Figures 1 and 2, and are denoted by 
Lx * (1 — N). One finds that the corrections due to these new contributions are very sizeable. In 
the x-range considered here the effect of the Lx-resummation is practically cancelled. Hence the 
calculation of also the less singular terms in the higher-order splitting functions is indispensable 
for arriving at sound conclusions on the small x evolution on the small- x polarized evolution. 



5 Conclusions 

We have investigated the effect of the resummation of terms of order c4 +1 In 2 ' x, derived from 
the infrared evolution equations in ||, on the small-x behaviour of polarized singlet parton 
distributions and the structure function gi in deep-inelastic scattering. The comparison with 
the corresponding contributions obtained in the same order by complete NLO calculations of 
the splitting function matrix || shows the equivalence of both approaches in this limit up order 
a 2 s . Since the coefficient functions up to two-loop order contain only terms less singular in In a; 
in the MS scheme, the contributions oc a^ln 4 x in the three-loop MS splitting functions P^ 2 \x) 
have been predicted on the basis of this resummation. 

As a general result for SU(N C ), the off-diagonal elements of the matrix P(x, cx s )x-*o are 
found to be proportional by the factor —TpN^/Cp- The supersymmetric relation holds for 
the leading small-x terms of the anomalous dimensions in J\f = 1 supersymmetric Yang-Mills 
theory. Starting with order I = 1 even more constraining relations are found in this case and all 
anomalous dimensions are related at a given order in a s . 

The numerical analysis shows that the all-order resummation of the terms 0(a[ +1 \n 21 x) 
leads to very large corrections at small x. Terms less singular in lnx, being not calculated yet, 
can be expected to contribute in a very significant way even at the smallest x-values considered, 
x ~ 10 -5 , as in the case of leading and next-to-leading order. Even a full compensation of the 
effect obtained resumming the most singular terms can not be excluded. Hence solid conclusions 
on the small-x evolution of polarized singlet parton densities and structure functions can only 
be drawn if these terms are calculated. 
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Figures 

Figure 1: The Q 2 evolution of the polarized singlet quark and gluon momentum distributions xAS 
and xAg as obtained from the GRSV standard distribution at Ql = 4GeV 2 . The results are 
shown for the NLO kernels (full), the leading small-a; resummed kernels (dashed), and a modification 
of the latter by possible less singular terms discussed in the text (dotted curves). 

Figure 2: The x and Q 2 behaviour of the polarized proton and neutron structure functions g[' n (x, Q 2 ) 
as obtained from the parton densities in the previous figure. The notations are the same as in Figure 1. 
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